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1. Introduction 

Much of our knowledge about hadronic structure in terms of quark and gluon degrees of free- 
dom has been obtained from high energy scattering experiments. However, as discussed in the talk 
of Vanderhaeghen [|l|], there are still many unresolved issues in hadronic physics that need to be 
addressed, from both an experimental and theoretical perspective. This is one of the main moti- 
vations of the 12 GeV Jefferson Lab upgrade which aims to [Q]: search for exotic mesons; study 
the role of hidden flavours in the nucleon; map out the spin and flavour dependence of the valence 
parton distribution functions; explore nuclear medium effecs; and measure the generalised parton 
distribution functions of the nucleon. It is imperative that these and other exciting experimental 
efforts, such as those at COMPASS/CERN and FAIR/GSI, are matched by modem lattice simu- 
lations which, thanks to recent innovative computer and algorithmic improvements [||], are now 
capable of reaching light quark masses (mj^ < 300 MeV) and large volumes (>3 fm) [Q]. 

In this talk I will report on progress made in the past year (for reviews of results reported in the 
previous two conferences, see earlier reviews by Orginos [Q] and Hagler |^]) in lattice calculations 
of many different aspects of hadronic physics such as the electromagnetic form factors of N,K, p, A 
and transitions in Sec. ^ moments of structure functions in Sec. ^ axial coupling constants of 
baryons in Sec. ^ moments of generalised parton distributions (Sec. |5|) and distribution amplitudes 
(Sec. ^, disconnected contributions in Sec. ^ polarisabilities in Sec. |8|, and finally in Sec. ^ I 
summarise the current status of these topics and point out unresolved issues and directions for the 
future. 

2. Electromagnetic Form Factors 

The study of the electromagnetic properties of hadrons provides important insights into the 
non-perturbative structure of QCD. The EM form factors reveal important information on the in- 
ternal structure of hadrons including their size, charge distribution and magnetisation. 

A lattice calculation of the ^^-dependence of hadronic electromagnetic form factors can not 
only allow for a comparison with experiment, but also help in the understanding of the asymptotic 
behaviour of these form factors, which is predicted from perturbative QCD. Such a lattice calcula- 
tion would also allow for the extraction of other phenomenologically interesting quantities such as 
charge radii and magnetic moments. For a recent review see [^. 

2.1 Nucleon Form Factors 

Phenomenological interest in the electromagnetic form factors of the proton has been revived 
by recent Jefferson Lab polarisation experiments [^] measuring the ratio of the proton electric to 
magnetic, ji.'^P^ G^f\q^) / Gm\q^) , and Pauli to Dirac, F2{q^) / Fi{q^) , form factors. Based on per- 
turbative QCD [^, the asymptotic scaling behaviour of these ratios should be independent of 
(for Ge/ Gm ) or scale as 1 / q^ (for F2/F1), however these experiments showed that G^. / Gm decreases 
almost linearly with increasing q-^, while F2/F1 scales as 1/y^. Additionally, fits of proton and 
neutron data using phenomenologically motivated ansaze provide for the possibility of a zero cross- 
ing in the isovector electric form factor, Gg, around ~ 4.5 (GeV/c)^ [p^. 
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The electromagnetic form factors of the neutron aie also receiving plenty of interest at the 
moment since we know that it has charge zero, but how is its internal charge distributed and does 



it have a positively or negatively charged core [|11|]? Lattice calculations can provide insights into 
this distribution since lattice simulations of three-point functions are performed at the quark level, 
and hence they have an advantage over experiment in that they can directly measure the individual 
quark contributions to the nucleon form factors. 

On the lattice, we determine the form factors F\{q^) and F2{q^) by calculating the following 
matrix element of the electromagnetic current 



{p',s'\f{q)\p,s) =u{p',s' 



2m I 



u{p, s) , 



(2.1) 



where u{p, s) is a Dirac spinor with momentum, p, and spin polarisation, s, q = p' — p is the 
momentum transfer, mn is the nucleon mass and is the electromagnetic current. The Dirac (Fi ) 
and Pauli (F2) form factors of the proton ai^e obtained by using 7^^^ = ^uy^u — ^dy^d, while for 
isovector form factors = uy^u — dy^d. It is common to rewrite the form factors F\ and F2 in 
terms of the electric and magnetic Sachs form factors, Ge=Fi+q^ / {IniNYFi and G„, = F\+F2. 



If one is using a conserved current, then (e.g. for the proton) F^'^^ (0) = GY' (0) = 1 gives the 



electric charge, while Gm ' (0) = /X = 1 + k:^'^) gives the magnetic moment, where F!f^ (0) = K 



(p), 



is the anomalous magnetic moment. From Eq. (2.1) with see that F2 always appears with a factor 
of q, so it is not possible to extract a value for F2 at q-^ = directly from our lattice simulations. 
Hence we are required to extrapolate the results we obtain at finite q^ to q^ = 0. Form factor radii, 
r; = ^J'(rf), are defined as the slope of the form factor at q^ = 0. 

In Fig. [I] we see results for the isovector Dirac radius from several different fermion actions. 
RBC/UKQCD presented an update from the A'y^ = 2 + 1 domain wall fermion run on 24^ x 64 
lattices with a^^ = 1.729 GeV [ [l2| ] (red circles), while LHPC updated their mixed action (DWF 



valence, asqtad sea) results at their lightest pion masses [ ]13| ] (green right triangles). Additionally, 
LHPC have started running on the 32^ x 64 DWF configurations with a^^ 2.4 GeV generated by 
the RBC/UKQCD collaborations, and preliminary results from these runs are shown by the black 



upside-down triangles []13[]. 

These latest results are compared with earlier quenched and Nf 



2 Wilson [14, 13] and DWF 



[16, 17| results. We observe agreement between the different lattice formulations, while any dis- 
crepancies are an indication for systematic uncertainties, such as finite volume effects, discretisa- 
tion errors, etc. The overall pattern is typical of lattice results for ri , i.e. the lattice results lie below 
experiment with little variation as a function of m|. Investigations using chiral perturbation theory 
predict that these radii should increase dramatically close to the chiral limit [|4, 18]. Current results 
indicate that in order to see such curvature, one needs to perform simulations at niji < 300 MeV. 

During the conference, we also saw a preliminary analysis from the European Twisted Mass 
Collaboration using Nf = 2 twisted mass fermions with pion masses down to m;t ~ 313 MeV at a 
single lattice spacing, a = 0.089(1) fm [19], and results are forthcoming. 

Finally, QCDSF have been studing the ^^-dependence of the individual quark contributions to 
the nucleon form factors. In Fig. ^ we see some results for the ratio of the d- to w-quark contribu- 
tions to the proton's Dirac radius. Here we clearly see that rf > r" for all simulated quark masses 
(the same behaviour is seen for r2), indicating that the d (M)-quarks are more broadly distributed 
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Figure 2: Ratio of Dirac radii for u and cZ-quarks 
from QCDSF. Dashed hne indicates physical nij^. 



than u (J)-quarks in the proton (neutron). Note that discoimected contributions were not considered 
in this study. 

2.2 Accessing small Q^: Partially twisted boundary conditions 

On a lattice of spatial size, L, momenta are discretised in units of In/L. Modifying the bound- 
ary conditions of the valence quarks Yixk + L) = e'^''\if{xk), {k = 1,2,3) allows one to tune 
the momenta continuously p + d/L. Momentum transfer in a matrix element between states with 
initial and final momenta, pj + di/L and pf + df/L, respectively, then reads = {pf — piY = 

{ [Ef{pf, df) -Eiipi, di)]^ - [{pf + df/L) - {pi + di/L)f] , where E{p, d) = ^m^ + [p + 6/LY. 

F2 is particularly interesting since it cannot be measured directly at = to obtain magnetic 
moments. Hence it needs to be extrapolatd from finite which can not only increase the error, but 
can also introduce a model dependence into the result. As can be seen in Fig. || from the QCDSF 
collaboration [p^], results obtained by using partially twisted bc's (open blue symbols) help to 
constrain the extrapolation to = 0. 

Twisted boundary conditions, however, introduce additional finite volume (FV) effects ~ 
e^'""^, which were shown to be small for the pion form factor in the Breit frame but can 



be substantial for iso vector nucleon form factors [g3|]. In p4[ | it was shown that when tbcs are ap- 
plied only to the active quarks attached to the current, the FV corrections depend on an unphysical 
and unknown parameter. They also found that the FV corrections are largest for the magnetic form 



factor with small twists. Indeed, the partially twisted be results at small q^ in Fig. H appear to be 



suppressed compared to the overall fit, which is the expectation from [ ]23| , |24[ ]. 
2.3 Large 

Lattice calculations suffer from noise at large Q^. This typically restricts the range of available 
momentum transfers to < 4GeV^. Earlier attempts by the LHP Collaboration to access the 
electromagnetic form factors out to 2^ ~ 6 GeV^ in the Breit frame, (/?' = —p = 2nn/L) p6| ] 
revealed that the relative error in F^{Q^), at fixed pion mass increases as n^. They found that in 
order to achieve a point at 2^ Rr: 6 GeV^ with a relative error of 30%, they would have to increase 
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Figure 3: Pauli form factor, F2{q^) together with an 
extrapolation to = 0. Open blue symbols indicate 
results using partially twisted bc's. 
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Figure 4: Fi{Q^) from [ p7[ obtained using a vari- 
ational analysis for nij^ ~ 480 (pink), 720 (purple), 
1 100 (grey) MeV. Dashed curve indicates a fit to ex- 
perimental data. 



the statistical accuracy by at least a factor of 50. Furthermore, to compound the difficulty, it was 
observed that the relative error in the iso vector Dirac form factor increased with ~ 

This has led the JLab group to attempt a study of these form factors using variational methods 
[27]. Their initial quenched study is performed on a 16^ x 64 anisotropic lattice = 3) using the 



Wilson gauge action and clover fermion action with 3 quark masses corresponding to pion masses 
of 1100, 720 and 480 MeV. 

To extract ( |2.1[ ) from a lattice 3-point function, accurate knowledge on the overlap factors and 
masses is required. Usually these are cancelled by constructing a ratio of 3pt and 2pt functions, 
however the drawback here is that often one needs to use a 2pt function with large momentum at 
large Euclidean times, which introduces additional statistical noise. An additional problem could 
arise if a smeared source is used that has been tuned at p = 0, but may not be ideal at large p. To 
circumvent these issues, the JLab group use three different choices of gaussian smearing and then 
solve a generalised eigenvalue problem to obtain the overlap factors and masses from the two point 
functions. These are then used in the 3pt correlator to solve for the form factors. 

Their preliminary result for Fi{Q^) using this method is shown in Fig. ^. Encouragingly, 
we see that they are able to find a clean signal up to ^ 5 GeV^. Simulations with dynamical 
fermions and lighter quark masses are now starting. 

2.4 Nucleon-/'n (Roper) Form Factors 

Using the methods outline in the previous section, the JLab group have performed a quenched 



study of the transition form factors of the ground-state nucleon to its Pi i excited state [ |28P 



{N2\j^{q)\Nl)=mAp') 



"M (P) 



(2.2) 



which are extracted from the correlators using a variational analysis. 

Results for (Pn \p) from this quenched initial study are shown in Fig. ||, where we clearly 
see that it is possible to obtain a signal in such channels. While the behaviour of the results is dif- 
ferent to that of the experimental data, this is probably a result of the heavy quark masses currently 
being used. Investigations are now under way with dynamical fermions and lighter quark masses. 
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2.5 Pion Form Factor 
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The pion form factor, Fji{Q^), has recently 
received a surge of interest from several lattice 
groups. This is partly due to the fact that the 
pion is the easiest hadron to study on the lattice, 
making it the perfect candidate for testing new 
techniques (e.g. twisted boundary conditions, all- 
to-all propagators). In addition to this, F-^{Q^) 
is an interesting quantity to study phenomeno- 
logically since its asymptotic {Q^ —>■ oo) normal- 
isation is known from n jJ. + V decay, and 
hence it allows us to study the transistion from 
the soft to hard regimes. At low the Fji{Q^) 
is measured directly by scattering high energy pi- 
ons from atomic elections [p^], however measure- 
ments at high require quasi-elastic scattering off virtual pions | |30| ] which leads to a model de- 
pendence in the extraction of the form factor from experimental data; a source of systematic error 
not present in a lattice calculation. 

Recently, RBC/UKQCD have used stochastic propagators with a single spin/colour source 



i 5 



2 3 
e'(GeV") 



Figure 5: Proton-fn factors from a quenched 
study [E8|. Masses as in Fig. 0. 



calculated using the so-called "one-end-trick" piQ, together with twisted boundary conditions [ |21| ] 
to calculate the Fji{Q^) at small values of the momentum transfer [32]. The results from this study 



are presented in Fig. ^ Here the smallest momentum transfer available on this lattice using periodic 
he's is denoted by the vertical dashed line and the results for F^{Q^) for these by filled circles. The 
results obtained using twisted boundary conditions are given by the triangles, and we clearly see 
that they smoothly fill the gap between the first fourier momentum and = 0. 

Using their results at the smallest 3 values of Q^, the authors computed the pion charge radius 



and, using the NLO expression from ChPT 03 3| ] 

12K 



('■;r)su(2) 



NLO 



/2 SK^f2 



mt 
log^ + 1 



(2.3) 



and the value of the pion decay constant in the chiral limit [34], are able to determine the LEC, 
Zg(mp) = —0.0093(10). Evaluting the expression using the physical m^c gives (r|) = 0.418(28) fm^, 
compared with (r^)exp = 0.452(11) fm^. 

The ETM Collaboration are also using the combination of stochastic propagators and twisted 
bc's to determine the pion form factor (and subsequently, (r|)) on their Nf = 2 twisted mass config- 
urations [35]. They have calculated (r^) at a number of pion masses and extrapolated their values 
to the physical pion mass using 2-loop ChPT [36], obtaining (r|) = 0.396(10) fm^. 



We also saw an updated analysis from the JLQCD collaboration with doubled statistics [p7|]. 
They are using all-to-all propagators to calculate Fji{Q^) with Nf = 2 overlap fermions on a 16^ x 
32 lattice with light quark masses down to a sixth of the strange quark mass. By comparing NLO 
and NNLO ChPt, they find the two-loop contribution to (r|) to be significant in the simulated 
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Figure 7: Comparison of latest lattice results of 
(r^) with the experimental value [ |3^ (dashed lines). 



region. Hence they performed a joint two-loop fit to (r|) and {r^g) (pion scalar form factor, see 
Sec. |L2l| ), from which they obtain (r|) = 0.404(22) (22) fm^. 

In Fig. 0, we present the current status of lattice determinations of (r|) by comparing the latest 

While there is a slight 



results with earlier determinations | |38| ] and the experimental value 
scatter, the general trend of the lattice results, even after attempts at including chiral logs, is to lie 
low compared to the experimental value. Whether this can be explained by finite volume effects, 
discretisation errors, or even the application of ChPT at such large quark masses, will require 
further investigation. 

2.6 A Electromagnetic Form Factors 

The matrix element of the electromagnetic current between spin-3/2 states has the form (c.f. 



Eq. (2.1) for spin- 1/2 states) 



(A(^',/)|/|A(^,.)) 



Eb{p')Eb{p) 



(2.4) 



where Ua{p,s) is a Rarita-Schwinger spin- vector, Mb is the mass of the decuplet baryon and 



aiiq')r + ^ip'^+pn 



4m| 



IniB 



■{p'^+p^ 



(2.5) 



The parameters ai, a2, c\ and C2 are independent covariant vertex function coefficients. For de- 
cuplet baryons, there are four multipole form factors, Geo, Ge2, Gmi, Gm3, which are defined 
in terms of ai, c?2, c\, ci, and are referred to as charge (£"0), electric-quadrupole [El), magnetic- 
dipole (Ml) and magnetic-octupole (M3) form factors, respectively. 

While the SO and Ml form factors give access to charge radii and magnetic moments in the 
same way as for spin- 1/2 baryons, the El and M3 moments accessible in spin-3/2 systems provide 
insights into the shape of decuplet baryons and have the potential to discriminate between various 
model descriptions of hadronic phenomena. 

The Adelaide group are in the process of finalising their analysis of the multipole form factors 



of the full baryon decuplet [40] in the quenched approximation. Their findings for the magnetic 
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Figure 8: Quenched Adelaide results for magnetic Figure 9: Gei form factors from the Cyprus group 
moments of the proton and A+. with exponential fits. 



moment of the A"*" is shown in Fig. |8] and the results are compared with earler proton results 
A simple quark model predicts that they should be equal, however, due to differing pion-loop 
contributions, the proton and A+ magnetic moments are expected to differ at the physical pion 
mass [^]. In fact, quenched ChPT predicts that the pion-loop contributions for the A+ come with 
an opposite sign to that of full QCD [43], and the results in Fig. |8| confirm this prediction. 

To confirm that this is a quenched artifact, it is important to perform simulations in full QCD 



at light enough quark masses. The Cyprus group have started to perform such simulations [ 144[ ] and 
they are currently simulating at pion masses down to ~ 350 MeV and so are now starting to enter 
the region where the quenched results "bend down". 

Both groups find that Ge2 is negative as shown in Fig. ^from the Cyprus group, indicating that 
A is oblate. The Adelaide group find that Gm3 deviates from zero only at small quark masses, while 
the Cyprus group only have a result at a single quark mass where they find that Gm3 consistent with 
zero, so it will be interesting to see if their results will deviate from zero as their results at smaller 
quark masses start becoming available. 

2.7 p Electromagnetic Form Factors 

The QCDSF collaboration has recently started an investigation into the electromagnetic form 



factors of the p meson [45]. The matrix element of the electromagnetic current between spin- 
1 states is decomposed in terms of three form factors, Ge{Q^), Gm{Q^), Gq{Q^). Of particular 
interest is the value of the quadrupole form factor at zero momentum transfer, Gq{Q^ = 0), which 
gives the quadrupole moment; a non-zero value would indicate spatial deformation of the p meson. 

Fig. [T^ shows the Sachs form factors for the smallest pion mass analysed {rriji f» 400 MeV). 
The electric form factor is fitted with a monopole ansatz and from the slope of the form factor at 
= 0, the charge radius is computed. After an extrapolation linear in to the physical pion 
mass, they find the preliminary result (rp) = 0.49(5) fm^, although it is reasonable to expect some 
chiral curvature to enhance this value. 

For the magnetic form factor, it is not possible to calculate directly at = 0, which is needed 
for the determination of the magnetic moment (^-factor). Hence, the results at ^0 are extrap- 
olated to 2^ = with a dipole ansatz. As discussed in Sec. 2.2, twisted bc's have the potential to 
help here. 
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Extrapolating linearly in m| to the physical 
pion mass gives gp = 1.6(1), smaller than the 
quenched result from the Adelaide group [ p^ ] and 



P = 5.29 K = .13620 Vol24^48 = 406MeV 



a study using the background field method [1470 . 

Similar to the magnetic form factor, the 
quadrupole form factor needs to be extrapolated 
from 2^ 7^ to 2^ = to obtain the quadrupole 
moment. The form factor is fitted linearly in Q^, 
although again twisted bc's will help to determine 
if this is a valid assumption, and the resulting mo- 
ments linearly in m|. The authors find a small 
negative result /i^ = —0.017(2) fm^, in agreement 



^ IJAf = 1.69 ± 0.25 


Gc(i), monopole fit. 

G'M(t), dipole fit 

linffiir fit 
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04£mf 


5 ^ ' S 


ii 
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;i 0.5 
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Figure 10: p form factors from a 24^ x 48 lattice 

with mn = 406 MeV, j3 = 5.29, ff = 0. 13620 
with 04 6|]. The negative result is interpreted as the 

p meson having an oblate shape, an interpretation enhanced by recent results using density-density 



correlators [|48|]. 

2.8 N ^ A Transition Form Factors 



We have seen in Sec. g^that there is emerging evidence that the quadrupole moment of the A is 
non-zero, indicating that the A is not spherically symmetric. The nucleon, on the other hand, being 
a spin- 1/2 particle, doesn't have a measurable quadrupole moment, however it still may possess 
an intrinsic quadrupole moment and thus also be spatially deformed. A possible way to search for 
such non-zero amplitudes is through the study of spin- 1/2 to spin-3/2 (yN —)■ A) transitions, which 
are also accessible in lattice simulations p9[ |50| , 

The matrix element for the vector N ^ A transition is defined in terms of three form fac- 
tors Gmi, Ge2, Gc2 which are known as the magnetic dipole, electric quadrupole and Coulomb 
quadrupole form factors, respectively. While the magnetic dipole is dominant, it is possible to 
search for non-zero quadrupole form factors by considering the following ratios measured in the 
lab frame of the A 



Rem{EMR) 



Rsm{CMR) 



(2.6) 



Gmi(G')' ' 2mA Gmi (G') 

Precise experimental data exists for these ratios and strongly suggest deformation of N and A. This 



has recently been confirmed in a full QCD simulation by the Cyprus group []51|], as we can clearly 
see in Fig. 11. While Rem = cannot be ruled out with the current precision on the hybrid run, 
RsM is clearly negative, in agreement with experiment. 



3. Moments of Structure Functions 

3.1 Nucleon Momentum Fraction, (x) 

Much of our knowledge about QCD and the structure of the nucleon has been derived from 
deep inelastic scattering experiments where cross sections are determined by its structure functions. 
Through the operator product expansion, the first moment of these structure functions are directly 
related to the momentum fractions carried by the quarks and gluons in the nucleon, {x)q_g, whose 
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sum must be Y,q{^)q + {^)g = 1- The scale and scheme dependence of {x)q and {x)g cancels out in 
the sum. 

Hence the quark momentum fractions are interesting phenomenologically and have been stud- 
ied on the lattice for some time. In fact, lattice studies of {x)q are notorious in that all lattice results 
to date at heavy quark masses exhibit an almost constant behaviour in quark mass towards the 
chiral limit and are almost a factor of two larger than phenomenologically accepted results, e.g. 
{x)^^^^ = 0.157(9), leading many a lattice practitioner to scratch their head and wonder "Will this 
thing ever bend down?", as predicted in [|^. 

To date, only connected contributions have been simulated to high precision, hence results 
are usually quoted for isovector quantities where disconnected contributions cancel. For the latest 
progress on disconnected calculations, see Sec. ^ 
Dynamical configurations are now becom- 
ing available at quark masses light enough to en- 
able calculations in the area where such bend- 
ing is predicted to set in. During the conference, 
RBC/UKQCD presented their findings from their 
Nf = 2 + 1 DWF configurations with pion masses 
as low as niji 
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330 MeV [12]. Results in the MS 



scheme at 2 GeV are shown in Fig. 12 and are 
compared with the latest results from the QCDSF 
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Figure 12: {x)„_d from RBC/UKQCD (DWF), 
QCDSF (Clover) and LHPC (Mixed) 



and LHP [54] collaborations. In this fig- 
ure we see excellent agreement between the older 
quenched [||, ||] and A^/ = 2 + 1 DWF runs and 
the Nf = 2 clover results, with the possible ex- 
ception of the lightest clover mass. This dis- 
crepancy may be attributed to a finite size effect 

(m^jL = 2.78), since these effects are expected to enhance (a) at light quark masses [^]. 

While we see agreement between the DWF and Clover results, we observe a gap between these 
results and those coming from the mixed action approach. Since the overall pion mass dependence 
is similar, this suggests that it is a renormalisation effect; a suggestion further enhanced when we 
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consider that the results from the mixed approach use (non-perturbatively improved) perturbative 



renormalisation []54[], while those from the other approaches use nonperturbative renormalisation 
of the operators involved. Of course, for this issue to be fully resolved, the mixed action results 
need to be renormalised nonperturbatively. 

The ;tQCD collaboration has also started an investigation of (a) using 16^ x32, Nf = 2 + \ 
Clover configurations from the CP-PACS/JLQCD collaborations with a = 0.1219 fm. Preliminary 



results from simulations at ~ 800 MeV were presented in [ ]58| , p9| ], with results from lighter 
quark masses forthcoming. 

RBC/UKQCD also presented results for the nucleon's helicity fraction, tensor charge and 



twist-3 matrix element, di [12|. 



3.2 Operator Product Expansion on the Lattice 

Moments of the nucleon structure functions can be expanded in the lattice regularisation as 

^(q^)=c^^\aq)A2(a)+c'^'^^-^A4{a) + ... higher twist , (3.1) 

r 

where q is the momentum transfer, a the lattice spacing, c^"^ the Wilson coefficients of twist, n, and 
A„ the reduced matrix elements. 

The leading twist matrix elements are nonperturbative quantities and can be studied on the 
lattice (an example of which we have just seen in the previous section). The corresponding Wilson 
coefficients, however, are usually calculated in continuum perturbation theory. Recently it has been 
shown that by applying the Operator Product Expansion to a product of electromagnetic currents 
between quark states, it is possible to determine the Wilson coefficients nonperturbatively [jSOj], 
allowing for a consistent treatment of the moments of structure functions. 

QCDSF are currently performing a quenched simulation on a 24^ x 48 lattice using over- 
lap fermions [^], which have the advantage that undesired operator mixings are suppressed by 
chiral symmetry and results are free of 0{a) artifacts. By considering two different momenta, 
q = ^(1,1,1,1) and ^(1,1,1,1), it was shown that preliminary results for Wilson coefficients of 
the 67 operators considered have the correct Bjorken scaling. Further improvements will involve 
using twisted bc's to access smaller momenta. With the full data available, a fully nonperturbative 
and consistent evaluation of the moments of nucleon structure functions will be possible. 



4. Baryon Axial Charges 

The axial coupling constant of the nucleon is important as it governs neutron j8 -decay and 
also provides a quantitative measure of spontaneous chiral symmetry breaking. It is also related to 
the first moment of the helicity dependent quark distribution functions, = Am — Ad. It has been 
studied theoretically as well as experimentally for many years and its value, gA = 1.2695(29), is 
known to very high accuracy. Hence it is an important quantity to study on the lattice, and since it 
is relatively clean to calculate (zero momentum, isovector), it serves as useful yardstick for lattice 
simulations of nucleon structure. 
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Figure 13: Scaling of gA with nij^L [t 
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Figure 14: gA, gzz, gzz from [|66p using a mixed 
action approach. 



4.1 gA 

The axial charge is defined as the value of the isovector axial form factor at zero momentum 
transfer and is determined by the forward matrix element 

{p,s\A"-'^\p,s) =2gAS^ , (4.1) 

where p is the nucleon momentum, and is a spin vector with = —mjj. 



gA has been studied in-depth for many years by the QCDSF [62] and LHP collaborations [63] 
and has been shown to suffer from large finite size effects. The RBC/UKQCD collaborations have 
recently calculated gA on their Nf = 2 + 1 DWF configurations [Q], where they observed the finite 



size effects to scale exponentially with mj^L [ p4( ] as seen in Fig. |13|for the DWF and Clover results. 



The ETM collaboration have also started simulations to measure gA on their Nf = 2 twisted 



mass lattices, and we saw a status report [|65|]. 



Finally, LHPC have a new simulated mixed action point at nij^ ~ 293 MeV and have also 



started to measure gA on the Nf = 2 + \ DWF configurations generated by RBC/UKQCD [|13|]. For 
the latter, measurements are being performed at three quark masses and two lattice spacings, but 
with similar volumes. Preliminary analysis indicates that results from the two approaches agree, 
indicating that effects due to unitarity violation in the mixed action approach is negligible. 

4.2 Axial Coupling Constants of Octet Baryons 

While there has been much work on the (experimentally well-known) nucleon axial coupling, 
there has been limited work on the axial coupling constants of the other octet baryons, which 
are relatively poorly known experimentally. These constants are important since at leading order 
of SU(3) heavy baryon ChPT, these coupling constants are linear combinations of the universal 
coupling constants D and F, which enter the chiral expansion of every baryonic quantity. 

Lin and Orginos []66|] have used DWF valence quarks on an Asqtad sea with mj^ ranging be- 



tween 350 and 750 MeV and their results for gA, gzi. and gzz are shown in Fig. 14 Fitting all 



three couphngs simultaneously using gA = D + F + ^„ Cj^ x", ^^s = ^ — + ^„ x", g^^ 
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F + Ln Cr^\ with X = (m| - ml)/{4nf-), they find 

gA = 1.18(4),tat(6)sys, gSS = 0.450(21 ).s,at (27)sys , gZZ = -0.277(15)stat(19)sys, (4.2) 

and D = 0.715(6)(29), F = 0.453(5)(19). Since there is little known from experiment for g^s and 
gzz, these results serve as a prediction and are in agreement with findings from ChPT and \aige-Nc . 

4.3 N* Axial Charges 

In the previous sections, we have seen results for axial couplings of ground state baryons. 
Recently, there has been an attempt to calculate the axial couplings of the two lowest lying, negative 
parity nucleon states, the A^*"" (1535) and A^*i (1650) [|^]. 

The authors have used the 16^^ x 32, Nf = 2 clover configurations from the CP-PACS collabo- 
ration with a = 0.1555(17) fm and nipjm,, = 0.804(1), 0.752(1), 0.690(1). In order to isolate the 
two negative parity states, they construct optimised source/sink operators from a combination of 
operators. In order to verify their method, they also calculate gA of the nucleon to compare with 
other determinations. They are able to see a signal and after extrapolating their results linearly in 
m\ to the physical pion mass, they find < 0.2, g^^ k. 0.55, which is consistent with the NR 
quark model. 

5. Generalised Parton Distributions 

Generalised Parton Distributions (GPDs) have received much attention, from both theory and 
experiment, in the past decade since they provide a solid framework in QCD to relate many different 
aspects of hadron physics, including form factors, parton distribution functions, impact parameter 
dependent PDFs and spin sum rules. The importance of these functions has led the QCDSF and 
LHP collaborations to perform lattice investigations of their moments [|^, where it is has 
been shown that the ^^-dependence of the generalised form factors associated with these moments 
flatten for increasing moment. This has the interpretation of a narrowing quark distribution in the 
transverse plane of a fast-moving nucleon as x^ ^ 1. 

Here we focus on the insights moments of GPDs provide into the spin structure of the nucleon. 

5.1 Spin Sum Rules 

It is now well known that quark spin carries only ~ 30% of the total spin of the nucleon, with 
the remaining ~ 70% coming from quark orbital angular momentum and glue. The total spin of 
the nucleon can be decomposed in terms of the quark and gluon angular momentum 

^=£/,(m') +7,(^2), (5.1) 

which is then further decomposed into the quark and gluon spin and orbital angular momentum 
contributions 

l = £iAI, + XL, + AG + L,, (5.2) 

^ 9 ^ q 

where AZ and AG are the standard gauge-invariant quark and gluon spin fractions, while the orbital 
angular momentum contributions are defined by Lq = Jq — Ar/2 and Lg=Jg— AG. The relation of 
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Figure 15: LHPC: AE'?/2 and U' (left) and evolution of L"^" with respect to the scale, (right) ni 



total angular momenta, J^ g, to the GPDs is due to Ji[]690 who showed that they can be expressed in 
terms of moments of GPDs 



J, 



0) 



where t = and A,B2o are matrix elements of the energy momentum tensor 



{P'\T^''\P) = U{P'] 



2m, 



20 



N 



(A^) + C2o{A'))U{P) 



niN 



(5.3) 



(5.4) 



Since A2jf(0) = (x)^'^ are simply the quark and gluon momentum fractions, we have by momentum 
conservation 1 = L^'^foC^) +'^2o(^)' hence we have a sum rule for the anomalous gravitomagnetic 
moments = ^5^20(0) +^2o(0)- Here we stress that although the sum is scale and scheme inde- 
pendent and is equal to zero, for the individual Bjq, this is not necessarily the case. 



Most of the work towards a determination of B|q has been done by the LHP [ 54 ] and QCDSF 



[JTOp collaborations. This year, we have seen an update from LHPC for their simulations using 
the mixed action approach (left plot. Fig. [l^) and some preliminary results of a study using the 
Nf = 2 + 1 DWF configurations from the RBC/UKQCD collaborations The results in the left 



of Fig. 15 indicate that the signs of the spin (AZ^) and orbital angular momentum (Lg) contributions 
are opposite for each quark flavour. The same behaviour has been observed by QCDSF [^]. The 
lattice result L"+'' ~ is in strong disagreement with relativistic quark models and has led LHPC to 
search for scale dependence in L*, as suggested by [^]. As seen in the right of Fig. 15, they find that 
j^u-d changes dramatically at small and in fact changes sign, which may help to reconcile the 
lattice and the quark model results, which generically are valid at a low hadronic scale /X ^ 1 GeV. 
Although as the authors point out, the one-loop evolution used here is not perhaps not quantitatively 
reliable below 1 GeV. 

A potential improvement in the determination of from the lattice is in the extrapolation that 
is required from the simulated points at 7^ to the required point = 0. As seen in Sec. |2.2| , this 
can be achieved through the use of twisted boundary conditions, which is currently being explored 



by the QCDSF collaboration |25]. This may become particularly important at light quark masses 
when the data becomes noisier, and hence the extrapolation is poorly constrained. 
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5.2 Spin Asymmetries 

In the past couple of years, lattice calculations of the moments of GPDs have provided faci- 
nating insights into how quarks are spatially distributed inside the nucleon [^] and pion [73|. Of 
particular interest is the strong correlation between the transverse spin and coordinate degrees of 



freedom [74], providing evidence for a sizeable Boer-Mulders function, hf{x,k^) [75]. 

Recently, there has been an attempt to determine on the lattice (moments of) the Transverse 
Momentum Dependent PDFs (TMDPDFs), e.g. f^{x,k±^), h^{x,k±), which are important in semi- 
inclusive DIS (SIDIS). In order to obtain information on the dependence of these functions on the 
transverse momentum, ^j^, of the quarks inside a hadron, it is necessary to consider matrix elements 
{P\q{i)r'^ ^(0)|P), where the quark fields are separated by a distance, i, and 'W is a. Wilson line 
(to infinity and back). Of course, this is not possible on the lattice so instead one is required 
to consider a path of finite total length i separating the quark and anti-quark in the operator, as 



illustrated in Fig. 3a of [76]. 

The matrix element is then obtained from from 



C3pt('^)^sink)-P)r) 0<T<Osi 
C2pt(fsink)-P) 



{P\q{£)r^q{0)\P)ocAi 



■P)- 



(5.5) 



Choosing F to be 7^ gives access to A2,A3, while 7^75 gives Ag, A7, Ag. The £^-dependence of 
these functions is fitted with a double Gaussian. Moments of the TMDPDFs are then obtained via 
a Fourier transform 



/r 



dxfi{x,ki 



,(l)lat 



and similarly for g\j which is obtained from Ay. 

Information on the correlation between the 
intrinsic quark transverse momentum and the 
transverse polarisation of the nucleon can then be 
obtained by considering the combination 



l(^(l)M 



■S±_ (l)lat., 
8 IT {^i 



(5.7) 



which is shown in Fig. [T^ for a longitudinally po- 
larised M-quark inside a nucleon that is transver- 
sly polarised in along the x-axis, S± = {Sjc,0). We 
clearly see that the distribution is distorted along 
the ;c-axis. 



{2k) 



L2A2(|£±|,0) , 



(5.6) 
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Figure 16: TMDPDF for a w-quark in a nucleon 
that is polarised in along the x-axis 



6. Distribution Amplitudes 

Distribution amplitudes (DAs) describe the momentum-fraction distribution of partons at zero 
transverse separation in a particular Fock state, with a fixed number of constituents. They are 
essential for the determination of the hard contributions to exclusive processes, but being universal 
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hadronic properties, are process independent. Hence, they are important for calculations of form 
factors at large Q^, B-decays, and can be related to the Bethe-Salpeter wave function. 

DAs are defined as non-local matrix elements on the light cone, e.g. the leading twist pion DA 

{o\d{-z)r^,r5[-z,zHz)\7i+{p)) = if^p^ J\^e-'^^^^ ^ =x-x. (6.1) 

Results for moments of the light pseudoscalar meson distribution amplitudes have been presented 
by QCDSF [^] and UKQCD/RBC @ in the last couple of years. Here we will focus on some 
recent results for vector mesons and the nucleon. 

6.1 Vector Mesons 

For spin-1 mesons, there are two DAs, 0^(<^), as opposed to a single DA for spin-0 

mesons. The lowest moments of 0" (<^) are obtained from the local matrix elements 

(0|^(0)7{p D^} s{0)\V{p,X)) = mvfvP{pe^^^{^')\ , (6.2) 
(0|^(0)7{p B^,Bv} qmVipM = -imvfvP{pPi,e^^^{^^)l , (6.3) 

where my and fy are the mass and decay constant, respectively, of the the vector meson, V , and 
£^ is a polarisation vector. The moments, {^")\ are extracted by constructing ratios of lattice 



two-point functions [ [79| , |80| | and the bare lattice results are then renormalised. 

In Fig. 17 we see some preliminary results from the RBC/UKQCD collaborations for (i^ ^)^* 
and (i§^)p calculated with Nf = 2 + \ DWF configurations with 4 values of the light quark mass 
and 2 volumes []80|]. The results indicate that there are no clear signs of finite volume effects. 
After renormalising perturbatively (although in [ [80| ] they also presented a status report on their 
nonperturbative renormalisation programme and the results should be finalised soon) to the MS 
scheme at jx^ = 4GeV^, they find 

(^)|, «0.0359(17)(22) (^^^J « 0.240(36)(12) « 0.252(17)(12) , (6.4) 



which compare well with the preliminary results from QCDSF [ [79[ ] {^)\* ~ 0.036(3), (i^)^. ~ 
0.030(2). These results show the SU (3) y-breaking effects in the K* DAs in a similar way to that 



observed for the K DAs in 1 77, 78]. 
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0.005 




Figure 18: Nucleon decay constant, fy, (top/left) and difference between ratios of moments (bottom/left) in- 
dicating an asymmetry between (j)^^ and 0"^*'. Right: Barycentric contour plot of the leading-twist nucleon 
distribution amplitude at = 2 GeV. 



6.2 Nucleon 

For the proton, there are three distribution amplitudes, V,A,T. In a similar way to the case of 
mesons above, their moments (v''"" , A''"" , T'"'") can be obtained from hadron-to-vacuum matrix 



elements of local operators |81]. It is useful to construct the combination, 0'™" = — a''™ + 



2T'"'"). In the asymptotic limit, q){xi,Q^ ^ oo) = 120x1X2X3 and we have 0^°" = 0°"' = 0°'" = i, 
0200 _ 0020 _ 0002 _ 0110 — 0101 _ 0011 _ hcncc it is useful to look for asymmetries, such 
as 0100 -0010. 



QCDSF have calculated first two moments [82] using an improved constrained analysis which 



considers ratios of correlators [82] together with nonperturbative renormalisation of the appropriate 



3-quark operators [^]. By considering the difference between two such ratio, as shown in Fig. [18|, 
the asymmetry is pronounced and increases as one approaches the chiral limit. 

These asymmetries are visualised in Fig. where the lattice moments have been used in a 
polynomial expansion of the full nucleon DA. Here xi 2,3 refer to momentum fractions of the three 
quarks in the proton and the asymmetries indicate that the M-quark with spin aligned with proton 
spin has largest momentum fraction (xi). Interestingly, the asymmetries are less pronounced than 



for QCD sum rules [84] and other phenomenological determinations [85]. 



7. Strange Quarks in the Proton 

The determination of the strange quark content of the nucleon offers a unique opportunity to 
obtain information on the role of hidden flavour in the structure of the nucleon. Since the nucleon 
has no net strangeness, the strangeness contribution to the total charge of the nucleon must be zero, 
i.e. (0) = 0. However, there is no such simple constraint on either the sign or magnitude of 
the strangeness contribution to the magnetic moment, G^(0). Additionally, the strangeness charge 
radius may also be non-zero. While the latest experimental results [ p6| ] suggest that the strange 
form factors of the proton are consistent with zero, forthcoming experiments at JLab and Mainz 
will further clarify this picture. Additionally, the strangeness contribution to the total spin of the 
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the contributions from the individual quark sectors ^-^ ' 
from a mixed action simulation. 



nucleon is poorly determined. Hence there is an opportunity for lattice simulations to make an 
important contribution to the current understanding of the role of strange quarks in the nucleon. 

7.1 Indirect strangeness 

An indirect method for determining the electromagnetic strangeness form factors has been 



proposed over the last couple of years by the Adelaide group [87]. By combining charge symme- 
try constraints with chiral extrapolation techniques, based on finite-range-regularisation [^], and 
low-mass quenched-QCD simulations of the individual quark contributions to the charge radii and 
magnetic moments of the nucleon octet, precise estimates of the proton's strange electric charge 
radius and magnetic moment were obtained. 



Recently, Lin & Orginos [ |27[ ] have followed this procedure using results from a mixed action 
simulation (DWF valence on Asqtad sea) with pion masses in the range m^t = 350 — 750 MeV. Their 
findings for the individual quark contributions to the charge radii of the nucleon octet indicate that 
the contribution from the heavier strange quark is smaller than those from the light quarks (Fig p^ , 
in agreement with quenched results pi]]. 

After taking these mixed action results and following the Adelaide method, Lin & Orginos 
findate^^O.lGeV^ 

G'm = -0.082(8)(25) , G| = -0.00044(1)(130) , (7.1) 



which is in excellent agreement with earlier findings [87] and recent JLab experiments [86]. 



7.2 Direct strangeness and other disconnected 

Direct lattice calculations of the strangeness content are computationally demanding and are 
renowned for suffering from large statistical noise. Recent advances in computing power combined 



with technical innovations, such as all-to-all propagators [89], have led to a renewed interest in 
direct determinations of disconnected quantities, such as strangeness in the nucleon. This year we 
have seen the progress being made in this area from several groups using a variety of different 
methods. 
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The ;tQCD collaboration 59] have started a simulation to determine the gluonic and 
strange quark momentum fractions of the nucleon, {x)g, (x).,, and the strangeness magnetic mo- 
ment using Nf = 2 + 1 dynamical clover configurations from the CP-PACS/JLQCD collaborations. 
They use Z(4) stochastic noise sources combined with an unbiased subtraction from the hopping 



parameter expansion (HPE) []91[], and multiple (up to 32) sources. By summing over the operator 
insertion times, they can then fit to the slope, as shown in Fig. Their preliminary findings sug- 
gest that the strange-to-Ught momentum fraction ratio, {x)s/ {^{{x)a + = 0.857(40), which is 
slightly larger than the CTEQ value, 0.27 < r < 0.67. Additionally, {x)g is studied using the over- 
lap operator to construct F^y [ ]92| ] in quenched QCD; since ultraviolet fluctuations are expected to 
be suppressed due to the exponentially local nature of the overlap operator. For {x)g, a signal is 
obtained with ~ 3a accuracy, however renormalisation is required. 

Another group that is making substantial progress is the Boston group [^3|]. They are using 
stochastic sources with maximum dilution with Nf = 2 Wilson fermions and mji f» 400 Me V. By 
using vacuum-subtracted currents, e.g. V — {V), they fit the three-point function directly using 
input from the two-point functions to calculate G^^{q^ = 0), G\{q^ = 0) = A^. They find that they 
are able to determine As with 30% errors and their result for is consistent with zero. From the 
scalar form factor, they obtain the result fys = "''^m"^^^ ~ 0.48(7)(3). 



There was an update [ ]94| ] on work outlined in [^] to calculate As and {N\ss\N) using various 
noise reduction techniques such as HPE, truncation solver method, truncated eigenmode approach 
and dilution so that the stochastic source is only defined on a single timeslice. They find a reduction 
in the stochastic variance at fixed cost of around 25-30. Disconnected loops are currently being 
calculated using Wilson propagators on a staggered sea for three valence quark masses and two sea 
quark masses, with plans to move on to a full Nf = 2 + 1 simulation in the near future. 

7.2.1 Scalar Form Factor 

An analysis of the chiral behaviour of the scalar radius of the pion, (r|), can lead to a deter- 
mination of the LEC £4 and it is expected to have an enhanced chiral logarithm as compared to 
the vector radius discussed in Sec. Hence it is a good place to search for chiral nonanalytic 
behaviour in the chiral regime. However, such a calculation would need to take into account of the 
disconnected contribution to the form factor, and as a result it has received little attention to date. 

As mentioned in Sec. the JLQCD collaboration are computing all-to-all propagators on 



their Nf = 2 overlap configurations [37]. These all-to-all propagators allow them to compute the 
scalar form factor of the pion, including the contributions coming from disconnected diagrams. 

From the slope of this form factor, they calculate the scalar radius of the pion, (r^)|, which is 
shown in Fig. ^ as a function of m^. Also shown in the plot is the result of a combined NNLO 
ChPT fit to (r^)^, {r^)^ and cy (see Sec. 2^), where we can clearly see the predicted chiral cur- 



vature at light quark masses. It will be interesting to see if this can be confirmed as results be- 
come available at masses below 300 MeV. After extrapolating to the physical pion mass, they find 
(r^)^ = 0.578(69) (46) fm^, in agreement with experiment. 

8. Background Field & Polarizabilities 

All of the results presented in the previous sections have been obtained using sequential source 



19 



Investigations ofhadron structure on the lattice 



James M. Zanotti 




methods, where a current with a particular momentum is inserted into one of the quark propaga- 
tors used to construct the hadron. Static quantities are then obtained from a fit to the momentum 
transfer, q^, dependence of the resulting form factors. An alternative method is to consider lattice 
simulations in the presence of constant electric, E, and magnetic, B, fields, and by studying the en- 
ergy shifts in a hadron as a function of the field strength, it is possible to extract not only magnetic 
moments, but also electric and magnetic polarisabilities, and pM- 

On a finite volume, discontinuities can occur as the quark crosses the boundary of the lattice 
unless the fields are quantised as qa^B = Inn/L. However, this means that for the volumes being 
used in current lattice simulations, the applied fields are so large that nonlinearities can arise (and 
dominate) in the B-dependence of the masses and possible distortions in the particles themselves. 



Aubin et al. IjO^] studied the finite volume effects of the magnetic moment of the A baryon. 
They perform a test with quenched lattices for two different spatial volumes ((1.6fm)^ and (2.4fm)^), 
they implement a "patch" to the field by adding thex-link modification A^(L— l^y^z^t) = —aBLyS^x 
ii X = L — \, resulting a field that is quantised in units of iTzjlJ. They find that when using un- 
patched data with qa'-B <2n/L, large finite size effects are seen on the small volume. However, 
after patching, so long as simulations are performed close to lizjlJ, reliable results are obtained, 
even on the small volume. 

Having justified their method, they then proceeded to use 2 + 1 flavours of stout-smeared 
clover fermions on anisotropic lattices with two spatial volumes and patched magnetic fields. One 
light quark propagator {m-,i ^ 366 MeV) and one strange quark propagator to obtain results for the 
magnetic moments of A++'+'^'° and Q. . Results for A^ and Q. are consistent with experiment, 
with the accuracy of Q. comparable to experiment. 

Tiburzi et al. [0] calculated for both neutral and charge hadrons using clover fermions 
on DWF sea as a test run, with valence DWF to follow. They also showed the benefit of patching, 
although here it is referred to as "including transverse links", by showing that it is possible to 
remove "spikes" in, e.g. the pion's effective mass, by including links such as mentioned above. The 
electric polarisability, «£, of hadrons are then determined by examining the quadratic dependence 
of the energy as a function of field strength, as in Fig. 22. 

For charged hadrons, one must also take into account the sum of the Born couplings to the 
particles total charge, leading to a modification of the time-dependence of the particles two-point 
function. As a result, for charged hadrons are extracted from the exponential time behaviour of 
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the two-point functions by fitting them with exp{—Et — (Q^S'^t^)/ (6M)). While errors obtained in 
this initial study are fairly large, a couple of particularly interesting results are the ratio a^^ / , 
which is found to be in good agreement with expectations from one-loop chiral perturbation theory, 
i.e. it scales as niK/mji, and «£■ of K*^ and K*^ which are found to be negative. For the proton and 
neutron, they find a| = 3.6(1.3) and a| = 8.8(5.9) xlQ-'^fml 

Alexandrou [98] showed that an earlier calculation of [99] can be improved by simulating 
with an exponential background electric field rather than the linear field used in []99|]. 

Finally, in [ ]100| ] we saw some preliminary results for j8^'" on Nf = 2 clover configurations 
from the CP-PACS collaboration with three lattice spacings but constant physical volume, and 
quark masses in the range 0.547 < mji/mp < 0.8. 



9. Conclusion & Outlook 



Due to recent computer and algorithmic improvements, lattice calculations of hadronic quanti- 
ties are now becoming available at pion masses as low as nijc ^ 250 MeV, and it is not unreasonable 
to expect that soon simulations will be performed close to the physical pion mass. However, as 
we have seen in, e.g. gA, finite size effects (FSE) are starting to become a serious issue. As a 
result, many groups are now planning future simulations on volumes as big as (4fm)^, in order to 
minimise these effects, although corrections from ChPT will still probably need to be taken into 
account. 

This year we have seen an impressive amount of progress in many different hadronic quanti- 
ties, providing fascinating insights into the structure of hadrons. From the slope of the electromag- 
netic form factors, charge radii are now being computed for hadrons such as tt, p, A'^, A in a region 
where we expect to see dramatic chiral curvature towards the physical point. However, as these 
radii are an indication of the size of a hadron, as mentioned earlier, FSE need to be considered 
carefully. 

The scaling of hadronic form factors is now receiving an increasing amount of attention. 
In particular, twisted boundary conditions are providing access to small Q^, but there also is work 
underway to attempt to probe the large Q region (> 4 GeV^). The small region is also an 
interesting place to study the Dirac and Sachs electric form factors of the neutron. The results that 
are now becoming available at small are not only able to help constrain static quantities such as 
charge radii and magnetic/quadrupole moments, but also the value of the generalised form factor 
Bioi^^), which at = provides the value of the anomalous gravitomagnetic moment, which is 
important in Ji's angular momentum sum rule. 

Lattice calculations of moments of generalised parton distributions are providing insights into 
the different quark contributions to the nucleon's spin and angular momentum, and current results 
indicate 7„ 46%, Jd ~ 0, L„+rf ?a 0. These moments are also providing evidence for non-trivial 
transverse spin densities in the pion and nucleon. 

Simulations with zero momentum transfer lead to moments of ordinaiy parton distribution 
functions, and include phenomenologically interesting quantities such as gA and (x). Here, FSE 
appear to playing an important role in the extraction of these quantities, especially for gA where we 
have seen FSE lowering the lattice results. While there appears to be a slight tension between the 
renormahsation of some of the lattice results for (x)„_j, the overall pattern seems to indicate that 
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we may now at last be entering the region where the results may start to "bend down" towards the 
phenomenological value. Although once again, FSE are predicted to become an issue close to the 
physical pion mass, so care will need to be taken to ensure this encouraging behaviour continues. 

Following the recent success of lattice calculations of the moments of the light pseudoscalar 
meson distribution amplitudes (DAs), there are now results becoming available for moments of 
vector meson and proton DAs. Results for the proton are providing evidence that asymmetries 
exist in the way the momentum of the nucleon is distributed amongst its constituent quarks, with 
the M-quark with its spin aligned to that of the proton carrying the most momentum. The results 
also indicate that the symmetries are less pronounced than in QCD sum-rules. 

While it is important to push these more "standard" hadronic measurements as far as we can 
with the new sets of dynamical configurations that are becoming available, it is also important to 
develop new ideas and techniques. This year we saw a number of innovative methods for accessing 
less well known quantities. 

By considering matrix elements of operators where the quark fields are spatially separated, 
moments of Transverse Momentum Dependent PDFs have been computed. From these moments, 
it has been seen that densities of longitudinally polarised quarks in a transversely polarised nucleon 
are deformed. 

To date, lattice calculations of hadronic quantities have neglected the contributions coming 
from disconnected diagrams, since these are notoriously difficult to compute. Recently, however, 
there has been a renewed interest in determining these disconnected contributions to investigate 
the strangeness and gluonic content of the nucleon and, in particular, their contributions to nucleon 
spin. Through the use of all-to-all propagators and various noise reduction techniques, it may now 
be possible to calculate some of these contributions with as small as 10% errors. 

Although background field methods have been around for a long time, they have only recently 
received a lot of interest, since traditionally the electromagnetic fields induced on currently sized 
lattices were too large. However recent developments show that it is now possible to consider 
fields that are a factor of L smaller As a result, magnetic moments and polarisabilities can now 
be extracted from these simulations with much more confidence. Additionally, it has recently been 
shown that it is possible to extract the electric polarisabiUties of charged hadrons from a lattice 
simulation, and preliminary results are promising. 

In summary, lattice simulations of hadronic observables have received a surge of interest over 
the past few years, such that we are now not only in a position to confirm experimental findings 
from a first-principles calculation, but also to provide predictions for, and in some cases to guide, 
future experimental programmes. 
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